Many-body orbital paramagnetism in doped graphene sheets 
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The orbital magnetic susceptibility (OMS) of a gas of noninteracting massless Dirac fermions is 
zero when the Fermi energy is away from the Dirac point. Making use of diagrammatic perturbation 
theory, we calculate exactly the OMS of massless Dirac fermions to first order in the Coulomb 
interaction demonstrating that it is finite and positive. Doped graphene sheets are thus unique 
systems in which the OMS is completely controlled by many-body effects. 



Introduction. — The diamagnetic properties of carbon 
allotropes such as diamond and graphite have always at- 
tracted a great deal of interest. Early on, it was shown 
experimentally that the diamagnetic susceptibility of 
graphite is large and strongly anisotropic [1 . These stud- 
ies stimulated an intense theoretical activity [2. The re- 
cent isolation of graphene (see Ref. 3 for reviews) has re- 
vitalized an interest in the magnetic properties of carbon- 
based materials. Graphene is a truly two-dimensional 
(2D) system composed of carbon atoms tightly packed 
in a honeycomb lattice. States near the Fermi energy 
of a graphene sheet are described by a massless Dirac 
equation which has chiral states in which the honeycomb- 
sublattice pseudospin is aligned either parallel to or op- 
posite to momentum [3] . McClure was the first one to re- 
alize that the Landau quantization for 2D massless Dirac 
fermions (MDFs) in a magnetic field perpendicular to the 
layer is very special (existence of a zero-energy Landau 
level) and to conjecture that this is the cause of the large 
diamagnetic susceptibility of graphite [2j [3] . 

More precisely, McClure showed that the orbital mag- 
netic susceptibility (OMS) of a noninteracting gas of 2D 
MDFs at a finite temperature T is given by [2j |4HZ] 



(o) 

^orb 



l 



24tt c 2 /c B T cosh 2 [/i/(2Jfe B T)] 



t^o g s gv e 2 v 2 



67T C 2 



S(e F ) , 



(1) 



where g s = g v = 2 are spin and valley degeneracies fac- 
tors, v is the Fermi velocity (which is independent on 
carrier density), c is the speed of light, \i is the chemical 
potential, and £p = fJ>(T = 0) is the Fermi energy. The 
zero-temperature OMS of graphene is thus infinite in the 
undoped limit, i.e. in the limit of zero carrier density. 
The second line of Eq. encodes another astonishing 
result [5j [6]. The OMS of graphene is exactly zero (at 
T = 0) if the Fermi energy is away from the Dirac point, 
i.e. if the system is doped. This situation seems to be 
really unique. 



Eq. ([I]) (and all the other studies [7] of the orbital prop- 
erties of MDFs we are aware of) heavily relies on a single- 
particle picture. The fundamental question we address in 
this Letter is the impact of electron-electron interactions 
on the OMS of a doped graphene sheet. Short-range re- 
pulsive interactions in a neutral Fermi gas [EJ [9] , for ex- 
ample, or Coulomb interactions in an ordinary parabolic- 
band electron gas [10] enhance the paramagnetic nature 
of the spin response (eventually driving the system to- 
ward a ferromagnetic instability [9j[10]) but are typically 
never strong enough to switch the sign of the orbital re- 
sponse from diamagnetic to paramagnetic. Orbital para- 
magnetism (OP), although possible in principle, is indeed 
a rare phenomenon in metals and semiconductors. For 
example, it has been shown that a 2D electron gas in 
a periodic potential exhibits OP when the Fermi level 
is sufficiently close to a saddle point of the band struc- 
ture [11 . Electrons in the proximity to a superconductor 
have been shown [12] to exhibit OP. Using diagrammatic 
perturbation theory up to first order in the Coulomb in- 
teraction we will demonstrate that the OMS of an inter- 
acting gas of 2D MDFs is finite and positive. Weakly- 
interacting doped graphene sheets thus represent unique 
systems with a paramagnetic orbital response of purely 
many-body origin. 

MDF model Hamiltonian and linear-response theory. — 
The (single-channel) Hamiltonian of a 2D gas of MDFs 
in the eigenstate representation is (Ti = 1) 
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where 6k,x = Xvk (A = ±) are band energies, S is the area 
of the system, p q = E fc ,A, M M x^k,q)c ] k _ q/ ^ x c k+q/2 ^ is 
the density operator, v q — 2ire 2 / (eq) is the 2D Fourier 
transform of the Coulomb potential, e is an average di- 
electric constant, and M\^(k^q) are matrix elements 
which can be found, for example, in Ref. [6j The many- 
body properties of doped graphene sheets depend [3j [13] 
on the dimensionless coupling constant (restoring ft for a 
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FIG. 1: Diagrams for the orbital magnetic susceptibility, a) 
The noninteracting bubble. Panel b) The first-order vertex 
correction. Panels c) and d) The two first-order self-energy 
diagrams. The wave- vector labeling of the propagators en- 
sures time-reversal symmetry and thus guarantees that the 
two self-energy diagrams give an identical contribution. 



moment) a ee = e 2 / (ehv), which can be tuned experimen- 
tally by changing the dielectric environment surrounding 
the graphene flake [14 j. The model ([2| requires the intro- 
duction of an ultraviolet cut-off k max [13] on the fe-sums. 
This must be done with great care, since, as discussed 
in Ref. 6 , the presence of fc max breaks gauge invariance, 
which must be restored in the calculations. 

The OMS of the system described by Xorb? can be 
calculated from [TO] 
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where x(o) is a short-hand notation for the static trans- 
verse current response function. The bare current- 
density operator j q for MDFs is proportional to the 
pseudospin-density operator [3]: j q = vfr q . In what fol- 
lows we will calculate x(o) U P to order q 2 in the long- 
wavelength q — » limit using many-body diagrammatic 
perturbation theory. 

Diagrammatic perturbation theory. — Fig. 1 shows the 
four diagrams that contribute to the transverse (q = qy) 
pseudospin response function x(o) to first order in the 



electron-electron interactions. Solid lines are noninter- 
acting Green's functions in the eigenstate representa- 
tion, G\(k,u) = [uj - e kyX + ^fc,A] -1 where r} k ,\ = 
f]X sgn(/c — &f,a) 5 with n = + . Dashed lines are electron- 
electron interactions. Here /cp,+ = &f = V / 47rn / '(g 8 g v ) 
is the Fermi wave number corresponding to an electron 
concentration n [15 and fcp ? _ = & max . 

In the static limit the zeroth-order bare-bubble dia- 
gram in panel a) gives x^°\q 0) = ao + &2<? 2 , with 
ao = — ^max/(47rv) and <22 = 0. As discussed in Ref. EJ 
the term O(q ) has to be subtracted away to restore 
gauge- invariance. The fact that a2 = originates from 
a perfect cancellation of intra- and inter-band contribu- 
tions. After the ad hoc regularization ao = 0, we find, in 
agreement with the second line of Eq. ([!]) , that the OMS 
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of the doped noninteracting system is zero: x 

Panels b), c), and d) in Fig. [I] show the remaining first- 
order diagrams: panel b) is the so-called "vertex correc- 
tion", while panels c) and d) contain two "self-energy" 
insertions. The evaluation of these diagrams up to or- 
der q 2 is lengthy and will be presented elsewhere. To 
achieve analytical progress it turns out to be particu- 
larly useful [16] to decompose the isotropic interaction 
Vk-k' in Eq. ([2| in angular momentum components as 
Vk-v = Em V m (k, k')e im ^-^ with <p k = kx and the 
pseudopotentials V m (k, k') defined in Eq. (10) of Ref. [TBI 
In what follows we will introduce dimensionless variables: 
wave vectors will all be measured in units of fcp> the pseu- 
dopotentials V m (k, h') in units of 27re 2 /(e/cF), and the re- 
sponse function x(q) m units of the single-channel MDF 
density-of-states ^(^f) = ^f/(27tv 2 ), £f = ^f- The ul- 
traviolet cut-off & max measured in units of fcp will be 
denoted by A = fc max /fcF- 

Modulo terms O (a ), which, as mentioned above, must 
be removed by hand to restore gauge invariance [6], the 
static transverse response function x(o) to first order in 
a ee reads 
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where the three dimensionless coefficients S n are given 
by the sum of vertex S^ ^ and self-energy Sn SE ^ contri- 
butions, E n si VC) + sl SE) : 
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It is easy to prove that f[ VC \k — > oo) — > 5/(128fc) + 
5(18C-7)/(576ttA; 2 )+0(A;- 3 ) and that f 2 VC) (k -> oo) -> 
-5/(128fc) + C(fc- 3 ). Here C ~ 0.916 is Catalan's con- 
stant. This implies that the sum f[ VG \k) + f 2 VG \k) 
decays like k~ 2 for large k: we can thus take the limit 
A — > oo in Eqs. (JsJ) and ^ finding a finite, cut-off inde- 
pendent result for the sum S^ -* + S^ '*. This property 
is identically shared by the sum + S 2 SE ' ) • Thus ver- 
tex and self-energy contributions are separately conver- 
gent in the ultraviolet limit. It is however crucial to take 
into account both contributions to have a finite result. 
Indeed, the quantities S^ -* and are separately 

strongly divergent: for the Coulomb potential we find 
/3 VC) (<? ->• 0) = -/ 3 (S E) (6> -> 0) = 1/(192tt6> 3 ) at small 
angles. The subleading terms in the Taylor expansions of 
f 3 VC \0) and f 3 SE \0) are also singular (oc 1/(9) but these 
singularities do not cancel out upon summing f 3 VG \0) 
with /f E) (6>) and are an artifact of first-order perturba- 
tion theory, which misses screening. These pathologies 
are commonly cured [10 by using a statically-screened 
Thomas-Fermi interaction: (restoring units for a mo- 
ment) v q = 2ne 2 /[e{q + ^tf)], 4tf = ^v^ee^F being 
the Thomas- Fermi screening wave number [T7] . 

Substituting Eq. Q into ([3| we finally find that the 
OMS is given by 
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where we have defined Af(a e 



-(S 1 + S 2 + S 3 )/(27r). 
Note that the final result (11) is formally beyond the 
first order in a ee since Thomas-Fermi screening intro- 
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FIG. 2: (Color online) The dimensionless quantity A/" = 
— (Hi + H 2 + H 3 )/(27r) as a function of the fine-structure con- 
stant a ee - Note that J\f(a ee ) > and that it depends weakly 
on a ee - Inset: a zoom for very small a ee - The horizontal axis 
is in logarithmic scale. The solid line represents the analytical 
weak-coupling result, Af(a ee — >• 0) — »• — ln(a ee )/(487r 2 ). 



duces non-linear dependencies on a ee . The quantity M 
is presented in Fig. [2] One can see, first of all, that 
the function J\f(a ee ) is positive, which corresponds to 
OP, and depends weakly on a ee . It is easy to show 
that Af(a ee -» 0) -» -ln(a ee )/(487r 2 ), implying that 
Xorb oc — a ee In (a ee ) in the weak-coupling limit [18 j. We 
remark that the paramagnetic nature of the orbital re- 
sponse, i.e. Xorb > 0, is stable with respect to changes in 
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the range of inter-particle interactions. Using the Hamil- 
tonian ^ with contact repulsive interactions of strength 
(i-e. v q = vo > 0), we find the analytical expression 

e 2 13 

* OTb = ^?256^°- (12) 
Discussion and conclusions. — Our findings can be tested 
experimentally in a variety of ways. The most direct 
one is to measure the thermodynamic magnetization as 
a function of an applied magnetic field. Ingenious setups 
for these type of measurements have already been suc- 
cessfully applied to various carbon structures [19] and 
to conventional 2D electron gases [20]. Spin and or- 
bital responses can be distinguished by applying a tilted 
magnetic field [21]. To enhance the magnetic response 
one can resort to graphene laminate [22] (actually, the 
first magnetic measurements in this system have already 
been done [23]). Another possibility is to use macroscop- 
ically large graphene films from carbon on copper foil, as 
demonstrated in Ref. [24] 

Interestingly, OP in doped graphene can be probed 
also by acoustic measurements. Indeed, mechanical de- 
formations of graphene are known to produce a pseudo- 
magnetic gauge field [25]. The induced vector potential 
A(r) is proportional to the deformation tensor Uij(r) and 
the pseudomagnetic field is given by the usual relation 
Bs(r) = V r x A(r). Using the basic equations of the 
theory of elasticity it is easy to prove that the coupling 
between deformations and electronic degrees of freedom, 
which occurs via A(r), leads to a renormalization of the 
shear modulus /i s — ^ ^ s (o) — Ms - ^IXorb^ 2 /^ 2 and thus of 
the transverse sound velocity uo 2 (q) = ^ s (q)q 2 /p, where 
Q2 is a coupling constant [26] and p is the mass density. 
The positive sign of Xorb thus implies a softening of the 
sound velocity with increasing q whereas diamagnetic re- 
sponse would result in the opposite behavior. 

In summary, we have shown that doped graphene 
sheets have a very intriguing orbital magnetic response. 
If electron-electron interactions are neglected, the OMS is 
identically zero. When electron-electron interactions are 
taken into account the OMS turns out to be finite. To 
the best of our knowledge, this is the first system we are 
aware of in which many-body effects control completely 
the orbital response. The sign of the OMS cannot be pre- 
dicted a priori. To first order in Coulomb interactions we 
have shown that it is positive. Weakly-interacting doped 
graphene sheets are thus many-body orbital paramag- 
nets. 
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